CHIRAL GLASS PHASE IN CERAMIC SUPERCONDUCTORS 
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A three-dimensional lattice of the Josephson junctions with a finite self-conductance is employed 
to model the ceramic superconductors. By Monte Carlo simulations it is shown that the chiral 
glass phase is stable in three dimensions even under the infiuence of screening. The nonlinear ac 
susceptibility and the compensation effect are also studied. The compensation effect is shown to 
be due to the existence of the chiral glass phase. In agreement with experiments, this effect is 
demonstrated to be present in the ceramic superconductors which show the paramagnetic Meissner 
effect. 
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I. INTRODUCTION 

Since the discovery ceramic high-Tc superconduc- 
tors it has been known that they may exhibit a 
glassy behavior reminiscent of the spin glass 0] . Re- 
cently, by the noise and ac susceptibility measure- 
ments Leylekian et. al. demonstrated that LSCO 
ceramics show the glassy behavior even in zero ex- 
ternal field . They also observed an intergranular 
cooperative phenomenon indicative of a glassy phase 
transition. This collective phenomenon may be in- 
terpreted in terms of the chiral glass (CG) picture 
using a three-dimensional lattice model of the 
Josephson junctions with a finite self-conductance. 
The order parameter of the CG phase is a "chiral- 
ity" which represents the direction of the local loo- 
supercurrent over grains The frustration es- 

sential to realize the CG phase arises due to random 
distribution of 0- and 7r-junctions with positive and 
negative Josephson couplings, respectively In 
this paper we review our recent results obtained by 
Monte Carlo simulations on the nature of ordering 
of the CG. One of the most important conclusions is 
that the screening effect does not destroy this phase 
in three dimensions (the vortex glass phase which 
may exist only in the non-zero external field is, 
in contrast, unstable under influence of the screening 
1^). It should be noted that more direct support of 
the CG has been reported by the ac susceptibility j|] 
and resistivity measurements on YBCO ceramic 
samples and by the aging effect in Bi2Sr2CaCu208 

Another issue of the present paper is to explain 
the so called compensation effect (CE) observed in 
some ceramic superconductors Overall, this ef- 
fect may be detected in the following way. The sam- 



ple is cooled in the external dc field down to a low 
temperature and then the field is switched off. At 
the fixed low T the second harmonics are monitored 
by applying the dc and ac fields to the sample. Due 
to the presence of non-zero spontaneous orbital mo- 
ments the remanent magnetization or, equivalently, 
the internal field appears in the cooling process. If 
the direction of the external dc field is identical to 
that during the field cooled (FC) procedure, the in- 
duced shielding currents will reduce the remanence. 
Consequently, the absolute value of the second har- 
monics 1x2 1 decreases until the signal of the second 
harmonics is minimized at a field Hdc = Hcom ■ Thus 
the CE is a phenomenon in which the external and 
internal fields are compensated and the second har- 
monics become zero. 

The key observation of Heinzel et al. is that 
the CE appears only in the samples which show the 
paramagnetic Meissner effect (PME) jl^Jl^ but not 
in those which do not. It should be noted that the 
intrinsic mechanism leading to the PME is still un- 
der debate ||l5|,|l6|. Sigrist and Rice argued that the 
PME in the high-Tj, superconductors is consistent 
with the d-wave superconductivity . On the other 
hand, the paramagnetic response has been seen even 
in the conventional Nb [p7|-|l9| and Al supercon- 
ductors. In order to explain the PME in terms of 
conventional superconductivity one can employ the 
idea of the flux compression inside of a sample. Such 
phenomenon becomes possible in the presence of the 
inhomogeneities or of the sample boundary . 

In this paper we explain the CE theoretically 
by Monte Carlo simulations. Our starting point is 
based on the possible existence of the CG phase in 
which the remanence necessary for observing the CE 
should occur in the cooling procedure. Such rema- 
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nence phenomenon is similar to what happens in 
spin glass. In fact, in the CG phase the frustra- 
tion due to existence of 0- and 7r-junctions leads to 
non-zero supercurrents. The internal field (or the 
remanent magnetization) induced by the supercur- 
rents in the cooling process from high temperatures 
to the CG phase may compensate the external dc 
field. 

Using the three-dimensional XY model of the 
Josephson network with finite self-inductance We 
show that in the FC regime the CE appears in the 
samples which show the PME. This finding agrees 
with the experimental data of Heinzel et al |T2[ . 



II. MODEL 

We neglect the charging effects of the grain and 
consider the following Hamiltonian 0,p| 



Jij cos(6'j - 0j - Ai-i) + 



2£ 



^ j2 , - ^ /' Mf)dr , (1) 



where 9i is the phase of the condensate of the grain 
at the i-th site of a simple cubic lattice, A is the fluc- 
tuating gauge potential at each link of the lattice, (f)o 
denotes the flux quantum, Jy denotes the Josephson 
coupHng between the i-th and j-th grains, C is the 
self-inductance of a loop (an elementary plaquette) , 
while the mutual inductance between different loops 
is neglected. The first sum is taken over all nearest- 
neighbor pairs and the second sum is taken over all 
elementary plaquettes on the lattice. Fluctuating 
variables to be summed over are the phase variables, 
9i, at each site and the gauge variables, Aij, at each 
link. <I>p is the total magnetic flux threading through 
the p-th plaquette, whereas <i>p^* is the flux due to an 
external magnetic field applied along the z-direction, 



^ext I if P is on the < xy > plane 

otherwise , 



(2) 



where S denotes the area of an elementary plaque- 
tte. In what follows we assume Jij to be an inde- 
pendent random variable taking the values J or — J 
with equal probability (± J or bimodal distribution), 
each representing and n junctions. 



III. EXISTENCE OF CG PHASE 

In this section we employ the finite size scaling 
technique to study the nature of ordering of the CG 
phase in three dimension. The external field is set 
to be equal to zero ($p^*=0). 

At each plaquette the local chirality is defined by 
the gauge-invariant quality. 



2^3/2 j^^ sin( 



A-ij), 



(3) 



where Jy ~ Jij / J and the sum runs over a directed 
contour along the sides of the plaquette p. The over- 
lap between two replicas of the chirality is 



— AT / , '^B ' 



(4) 



where Np is the total number of plaquettes. In term 
of this chiral overlap, the Binder ratio of the chirality 
is defined as follows 



gcG 



[< qt >]/[< ql >f). 



(5) 



Here < ... > and [...] represent the thermal and 
the configuration average, respectively, gcc is nor- 
malized so that it tends to zero above the chiral- 
glass transition temperature, Tcg, and tends to 
unity below Tea provided the ground state is non- 
degenerate. 
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FIG. 1. The temperature and size dependence of the 
Binder ratio gcc for C = 1. Inset is a magnified view 
around the transition temperature Tea ~ 0.286 ± 0.01. 
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Monte Carlo simulation is performed according to 
the replica exchange method |2^. Equilibration is 
checked by monitoring the stability of the results 
against at least three-times longer runs for a sub- 
set of samples. The free boundary conditions are 
employed. Fig. 1 shows the results for the system 
sizes L — 3, 4, 6, 8 and 10, and for the dimensionless 
inductance C — \ {C — {2tt / (l)o)'^ J C) . The num- 
ber of samples we used is running from 1500 to 100 
depending on the system size L. Obviously, all of 
the curves of gcG for L = 3,4,6 and 8 cross at al- 
most the same temperature, strongly suggesting the 
occurence of a finite-temperature CG transition at 
TcG = 0.286 ± 0.01. A more careful analysis |] 
shows that the CG phase exists if C is smaller than 
a critical value Cc, where 5 < £c < 7. Thus, the CG 
phase is stable against the screening if the latter is 
not strong enough. 



IV. COMPENSATION EFFECT 
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In order to study the CE one has to apply the 
external field H which includes the dc and ac parts 



H 



Hdc + Hac COs{LUt) 



(6) 



It should be noted that the dc field is necessary to 
generate even harmonics. The ac linear suscepti- 
bilty of model (1) has been studied [|j by Monte 
Carlo simulations. Here we go beyond our previous 
calculations of the linear ac susceptibility We 
study the dependence of the second harmonics as a 
function of the dc field. In this way, we can make a 
direct comparison with the CE observed in the ex- 
periments p2[ . 

The dimensionless magnetization along the z-axis 
mormalized per plaquette, m, is given by 



1 



Np(j)o 



E 

p£<xy> 



(7) 



where the sum is taken over all Np plaquettes on the 
< xy > plane of the lattice. The real and imaginary 
parts of the ac second order susceptibility x'2i^) 
x'-ii^) calculated as 



nh„ 



rh{t) cos{2Ljjt)d{Lut) , 



m{t) sm{2Lut)d{ujt) , (8) 



where t denotes the Monte Carlo time. The dimen- 
sionless ac field hac, dc field hdc are defined as follows 



FIG. 2. The second harmonics obtained after field 
cooling in a dc field hdc = 1 from T = 0.7 to T = 0.1. 
The temperature is reduced in steps of 0.05. At the 
lowest T = 0.1 the dc field used in cooling is switched oH' 
and the second harmonics are generated by applying the 
combined field (6) . The dc field is stepwise reduced from 
hdC = 1 to hdC = 0. The inductance is chosen to be equal 
to £ = 4. The arrows indicate the sense of the changes 
in the dc field. The results are averaged over 40 samples 
and are qualitatively the same as those presented in Fig. 
1 of Ref. [12]. 

For model (1) the PME appears clearly for hdc < 1 
So the largest hdc we take is 1. On the other 
hand, as mentioned above, the CG phase is found 
to exist below a critical value of the inductance Cc 
where 5 < £c < 7. One has to choose, therefore, an 
C which is smaller than its critical value. To be sure 
that we are in the CG phase we choose C = 4 and 
T = 0.1. 

The second harmonics have been obtained by em- 
ploying Monte Carlo simulations based on the stan- 
dard Metropolis updating technique. While Monte 
Carlo simulations involve no real dynamics, one can 
still expect that they give useful information on the 
long-time behavior of the system. We take L = 8 
and oj = 0.001. The sample average is taken over 20- 
40 independent bond realizations. We set hac = 0.1, 
corresponding to « 0.016 fiux quantum per plaque- 
tte. Smaller value of hac turned out to leave the 
results almost unchanged. 

Fig. 2 shows the dependence of second harmon- 
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ics |X2|, |X2| = yOtU^TMF, on hd, in the FC 
regime. Our calculations follow exactly the exper- 
imental procedure of Heinzel et al First the 
system is cooled in the dc field hdc = 1 from T = 0.7 
down to T — 0.1 which is below the paramagnet- 
chiral glass transition temperature Tcg ~ 0.17 
The temperature step is chosen to be equal to 0.05. 
At each temperature, the system is evolved through 
2x10'* Monte Carlo steps. When the lowest tem- 
perature is reached the dc field used in cooling is 
switched off and we apply the combined field given 
by Eq. (6). Using Eq. (7) and (8) we monitor the 
second harmonics (the technique for obtaing the ac 
susceptibility may be found in Ref . ) reducing the 
dc field from hdc = 1 to zero stepwise by an amount 
of A/irfc — 0.05. 1x2 1 reaches minimum at the com- 
pensation field hcom — 0.7 ± 0.05. At this point, 
similar to the experimental findings p^ , the inter- 
section of X2 a-nd X2 observed. This fact indicates 
that at Hcom the system is really in the compensated 
state. Furthermore, in accord with the experiments, 
at the compensation point the real and imaginary 
parts should change their sign jl^ . Our results show 
that x'2 changes its sign roughly at hdc = hcom- A 
similar behavior is also displayed by X2 but it is 
harder to observe due to a smaller amplitude of X2- 

In conclusion we have shown that the finite- 
temperature CG phase is not spoiled by the screen- 
ing. The CE may be explained, at least qualitatively, 
by using the CG picture of the ceramic superconduc- 
tors. The CE is shown to appear in the CG phase 
in which the PME is present but not in the samples 
without the PME 
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